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ABSTRACT 



We study preheating in a general class of supersymmetric hybrid inflation model. Supersymmetry leads to 
only one coupling constant in the potential and thus only one natural frequency of oscillation for the homo- 
geneous fields, whose classical evolution consequently differs from that of a general (non-supersymmetric) 
hybrid model. We emphasise the importance of mixing effects in these models which can significantly change 
the rate of production of particles. We perform a general study of the rate of production of the particles 
associated with the homogeneous fields, and show how preheating is efficient in producing these quanta. 
Preheating of other particle species will be model dependent, and in order to investigate this we consider 
a realistic working model of supersymmetric hybrid inflation which solves the strong-CP problem via an 
approximate Peccei-Quinn symmetry, which was proposed by us previously. We study axion production in 
this model and show that properly taking into account the mixing between the fields suppresses the axion 
production, yet enhances the production of other particles. Finally we demonstrate the importance of back- 
reaction effects in this model which have the effect of shutting off axion production, leaving the axion safely 
within experimental bounds. 
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1 Introduction 



Reheating in the post-inflationary Universe is a very important process, since it describes 
the formation of all matter and energy in the Universe today However, it has only 
relatively recently been realised that in certain cases the process of reheating involves a 
stage of explosive production of particles by parametric resonance which cannot be described 
by the usual methods of the elementary reheating theory 0. This stage of 'preheating' is 
a non-perturbative and out of equilibrium process and as such is a very involved and little 
understood proceeding. In fact there have been complex analytical methods of preheating 
developed for only the very simplest of single- field inflationary potentials @]. It is generally 
numerically easy to obtain results for the first stages of particle production, yet a complete, 
self-consistent analysis should include all effects of backreaction of produced particles and 
their rescattering. At present lattice simulations (such as in 0, §[) must be performed to 
fully study these effects. Nevertheless, much progress is being made in developing methods 
to cope with the rich topic of preheating 

For a while now, the advantages of using more than one scalar field to implement inflation 
have been well appreciated. Hybrid inflation models ^ are probably the most popular of all 
present inflationary models. They succeed in removing the need for fine-tuning of the inflaton 
potential by providing a natural means to end inflation. Most importantly however, they 
allow inflation to occur at a relatively low energy density using small scale fields avoiding 
strong CMB constraints, a task which is virtually impossible to achieve using only one field. 
Thus, hybrid-type models are attractive from the point of view of particle physics and in 
particular are naturally well suited to realisations in supersymmetry (SUSY) where there is 
an abundance of scalar fields and flat directions. 

Preheating in hybrid inflation has been studied PJ using the standard toy hybrid potential of 
two fields V{(f)i, 02), coupled to an extra, massless scalar field, x by the interaction ^{hfcpl + 
hl(f)l)x^. It is found § that over most of the parameter space, preheating is rather inefficient 
in hybrid inflation with little production of (pi, 02 particles although production of x particles 
is possible for certain values of hi,h2. However, the results concerning the production of (pi 
and 02 particles may change in a supersymmetric version of the model, where the potential 
is derived from a supersymmetric superpotential of the kind W ~ n4>i4>2 + ■ ■ such that the 
same coupling will give the quartic self-interaction for the 02 held and the interaction term 
between 02 and the inflaton 0i. Therefore, supersymmetry will ensure that the masses of 
both fields at the global minimum of the potential are not only the same order of magnitude 
but equal. This in turn will give rise to a non-chaotic evolution of the classical homogeneous 
fields at the beginning of the preheating/reheating era. The classical fields will begin to 
oscillate with relatively large, little damped amplitudes, which will favour production of 
particles occurring. 

In addition, in hybrid models inflation ends with a phase transition when the inflaton field 
reaches a critical value. At this point, the second field 02 is sited at a maximum of the 
potential with vanishing mass. Once the inflaton passes this point, the mass of 02 becomes 
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negative allowing both fields to roll down the potential towards the global minimum. It is 
clear that across the region of the spinodal instability (negative curvature) , production of 
(1)2 particles can be quite intensive. Therefore, previous to any oscillation of the homogeneous 
fields, preheating will begin with an initial burst of 02 production. If the amplitude of the 
oscillations is large enough, this will also happen in the successive oscillations every time 
the classical fields enter the region of the spinodal instability. Although this only directly 
affects the rate of production of 02 particles, the effect is propagated into the production of 
01 particles through the mixing term in the mass squared matrix. 

The layout of the paper is as follows: section 2 briefly describes the theory of preheating 
as applied to many scalar fields and we point out that the coupling between fields should 
in general be taken into account; in section 3 we study a general model of supersymmetric 
hybrid inflation, focusing on the issue of production of 0i and 02 particles; in order to 
motivate production of other particles, which is model dependent, in section 4 we review 
a particular model of inflation where the strong-CP problem is solved via a Peccei-Quinn 
symmetry which leads to the production of axions. It will be important to calculate the 
effects of preheating in this model since the number of axions produced during reheating 
must be tightly constrained in order to be cosmologically acceptable. The results of the 
numerical study of the first stages of preheating in this model are presented, and compared 
to those results without including mixing terms. The first effects of backreaction of produced 
particles is studied and the total number of each type of particle calculated. 



2 Theory of preheating with multiple fields 

We begin with a model of the early Universe described by a potential, \^(0a), in an expanding 
Universe with scale factor a{t). Amongst the fields 0q, in the potential, there are degrees of 
freedom corresponding to scalar fields in which all energy is concentrated during inflation. In 
the usually considered single-field inflation this is the slowly rolling inflaton, 0, with effective 
potential V{(j)). In the hybrid scenario there are two such fields important to inflation, one in 
which the potential is relatively flat during inflation which acts as the inflaton field, 0i, and 
one which sits in a false vacuum, 02, and causes the end of inflation in the phase transition 
to its true vacuum. All other fields are not important to inflation serving only perhaps to 
provide the constant vacuum energy during inflation. 

In order to investigate the reheating of the Universe and determine the extent of the effects 
of preheating, we must first study the evolution of these fields as inflation ends and beyond 
into the regime of coherent oscillations. 

We may write the full quantum fields as: 

(/)„(x,t) = 0„(t) + 50„(x,t), (1) 

where 50q(x, t) represents the quantum fluctuations about the homogeneous zero modes, 
0Q,(t), and necessarily satisfies (50a(x, t)) = 0. 
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Neglecting ensuing effects of particle production, we proceed by considering the motion 
classically, solving the equations of motion for the classical homogeneous fields 4>ait), but 
using quantum fluctuations ~ Ho/2it (where Hq is the Hubble constant during inflation) in 
the initial conditions at the end of inflation to initiate the spontaneous symmetry breaking. 

The equations of motion governing the behaviour of the background fields, 0a (t), are: 

0a + 3i/0a + dV/d(j)c. = 0, (2) 

where the Hubble constant is given by the Friedmann equation. 
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(3) 



and H = a/a. 



Here, the subscript a runs over all the classical scalar fields in the theory relevant to infla- 
tion and V = V{(j)a) is the effective potential of these fields. During inflation all energy is 
contained in this potential and the Universe is in a vacuum-like state with vanishing tem- 
perature, entropy and particle number densities. The inflaton field, 0, is slowly rolling with 
the friction term dominating in Eq. (|]) so that: 



-1 dV 



3H. 



(4) 



This also describes the relevant situation for each classical field at the beginning of the 
reheating era. 



The equations of motion for the quantum fluctuations of the scalar fields in linear perturba- 
tion theory are given by: 



0, 



(5) 



where S(j)ka{t) are the time dependent Fourier modes (in comoving momentum space k) of 
the quantum fluctuations 50a, 



(x,t) 



y2(27r)3/2 



(6) 



and al and cik are creation and annihilation operators. Thus, Eq. (^ describes the interac- 
tions between the classical fields, 0a, and the quantum fluctuations, 50a. When neglecting 
backreaction effects due to the particles produced, the Hubble constant, H{t), and the mass 
squared matrix, A4'^^{t) = d'^V/d(l)ad(f>f3, are functions of the classical fields only, as given by 
the classical equations of motion (§) and (^. It should be clear that the equations of motion 
for the quantum fluctuations will in general be coupled through mixing terms in the mass 
matrix A^^^ 7^ (for a ^ jS); in other words the rate of production of each separate particle 
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is dependent on the rate of production of the others, even at the initial stage of preheating 
when backreaction effects are not yet important. 



Changing to comoving fields, ipka = o?^'^5(t)ka^ these equations can be simplified further to: 

^fca i^ka + Mlsi^ks = 0, (7) 



a2 



where one usually neglects the pressure term, e = '^{H'^ +2a / a)il)k^ which is a small correction 
that approaches zero as the Universe approaches matter domination (a ~ ^2/3). Typically the 
transition between the vacuum dominated state of inflation and that of matter domination 
during reheating is fairly rapid, so this term will be negligible. 

The solutions of (^, ipkait)^ complex functions whose initial values at the end of inflation 
at time to are taken to be |p!0|: 



where, 

u.l = e/a' + ML. (9) 

The comoving occupation number of particles of wave number k is an adiabatic invariant 
defined by: 

= — + I^.J j-2- (10) 

An instability in the growth of a given mode, il^kay corresponds to an exponential increase 
in the occupation number of particles, oc exp (2/ifc„mt). rh denotes the characteristic 
frequency of oscillations of the classical fields, given by their mass at the global minimum. 
The efficiency of preheating of a given particle with a given mode k is measured by the 
growth parameter fik^. 

Of course to be fully self-consistent, the potential and its derivatives in the equations for 
the classical fields and Hubble constant, (|^), (^, should contain the full quantum fields in 
the theory, V{cf)a), taking into account the extra particles being produced. In the Hartree 
approximation this is done via (50^) 7^ 0, where the expectation value contains all modes: 

m = J ^-^\s<pkj'. (11) 

The backreaction on the classical fields will not be important in the first stage of the evolution 
when (5(/>^) — 0, but the modes will grow and eventually dominate over the other terms. 
The classical evolution will be altered and in turn this will affect the quantum equations (^) 
which also must include further terms corresponding to produced particles. 

As a first step we may neglect production of particles altogether and study the evolution 
classically. Once the behaviour of the classical fields in a specific model is given, we may 
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try to rewrite the evolution equation for each quantum fluctuation, Eq. (|^), in the form of 
a Mathieu equation 



i^ka + (Mk) -'2qacos2z)ijk^ = 0, (12) 

where dash denotes differentiation with respect to z, and z oc rht. The {A, q) plane is 
divided into instability/stability regions, such that in the unstable regions the solution grows 
exponentially, ipk oc e^''^, and this can be interpreted as production of particles for the 
associated mode. Therefore, once the parameters Aa{k) and qa are known for a particular 
model we can get some insight about which particles are more likely to be produced, without 
the need of numerical integration. 

In the next section we will apply the above receipt to a general super symmetric hybrid 
model. First we analyse the motion of the classical fields when they begin to oscillate, 
focusing afterwards on the implications this will have for the production of their own quantum 
fluctuations during preheating. 



3 Supersymmetric hybrid inflation. 

In the hybrid scenario with just two real fields, the potential is given by, 

V{<P, N) = ^ (iV^ - 2Af + ^^^AT^ + imj^^ , (13) 

where (f) will play the role of the infiaton. During inflation, the field N is trapped in its false 
vacuum, = 0, whilst the (p field slowly rolls along the potential from above the critical 
value 0c = iV2X/g)A towards the origin. Once reaches the critical value, the mass 
changes sign and becomes negative, allowing to roll away from its zero value towards the 
global minimum, which in the positive quadrant is given by, 

00 = 0, No = V2A. (14) 

The fields then oscillate about their true vacuum values ending inflation and beginning the 
preheating/ reheating phase. 

Integrating the classical equations of motion, we obtain the trajectory in the cf) — N plane. 
In many cases after the initial symmetry breaking, one field may oscillate with an amplitude 
much larger than that of the other field which is essentially fixed at its VEV. In this case 
almost all the energy is contained in just one field and reheating is then similar to the single 
field case. This will greatly simplify the investigation of preheating, although in general it 
will not be possible. Alternatively both fields may oscillate, and the motion may be chaotic. 
Having one or another situation will depend on the ratio of the coupling constants, A/t?^, and 
they have been analysed in P], the chaotic motion being typical of a potential with A ~ g^. 
However, a different scenario will arise in the supersymmetric hybrid models of inflation for 
which the exact relation A = (7^/2 holds. 
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In supersymmetric hybrid inflation, the potential in Eq. (13) can be derived from the 
superpotential, 

ly = k$(n2 - a2) , (15) 

where $ and N are complex superfields, and N in the potential being the real components 
of the associated canonically normalised complex fields, and k"^ = X = g'^/2. In addition, 
the inflaton mass will be given either by a soft SUSY-breaking mass or else generated 
by radiative corrections to the inflaton potential. In either case, it has to be smaller than 
the combination kA to ensure a potential flat enough during inflation to allow for the slow 
rolling scenario [12|, that is. 
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(16) 
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where prime denotes derivative with respect to 0. Moreover, the model should be able to 
reproduce the correct level of density perturbation, responsible for the large scale structure 
in the Universe, accordingly to the COBE anisotropy measurements. The spectrum of the 
density perturbations is given by the quantity 

32 VI 
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where e^, is defined roughly 60 e-folds before the end of inflation, and the COBE value is 



60r} 



Sh = 1.95 X 10"^ [|l|]. Assuming \r]\ < 1/60, then (p^ = cp^e' 
constraint, 

kA ~ 2.36 X 10^° {nm^f^^ GeV , 

or similarly. 



and Eq. ( |18|) yields the 
(19) 



kA ~ 1.27 X 10'^\r]\ . (20) 
The value |?7| ~ 0.01 is also a reasonable assumption in order not to generate too large a 



tilt in the spectral index, n = 1 + 2// — 6e, whose present observational constraint is 
In - 11 < 0.2. 



The superpotential in Eq. ([T5|) provides not only the couplings of the inflaton and hybrid field 
N neccessary to achieve the hybrid mechanism, but also a non-vanishing constant vacuum 
energy V{0) through the scale A, which is also related to the critical value 0c- Hybrid 
inflation models derived from this kind of superpotential are called F-term hybrid inflation. 
A different version of SUSY hybrid inflation is provided by D-term inflation [|15|, where the 
critical value and the constant potential energy are related to the Fayet-Illiopoulus term 
coming from an anomalous U{1) symmetry. Alternatively the constant potential ^(0) may 
owe its origin to the hidden sector of a supergravity (SUGRA) theory |jl6|. Nevertheless, 
the ultimate origin of ^(0) and the scale 0c is not relevant for the discussion of preheating. 
We only demand a common coupling giving the quartic interaction term for N and the 
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interaction between and A^, and the condition that the potential vanishes at the global 
minimum. Therefore, the SUSY hybrid potential can be written as: 

N) = \^\N' - 2<Plf + K^N^' + \ml<p\ (21) 

where the origin of (pc is left unspecified. The slow- rolling and COBE constraints, Eq.s 
( p!6l - pO| ), remain the same with the replacement 0^ = A. 



3.1 Classical field evolution 



Whatever the values of the parameters k, (pc and which satisfy the above constraints, the 
behaviour of the classical fields when they begin to oscillate will follow a common pattern, 
dictated by the form of the potential, and we can make some general observations which are 
independent of the particular values of the parameters chosen: 

i) The maximum attainable amplitude of the oscillations of the field, 0c = No/\/2, is 
comparable to (more precisely is a factor of \/2 smaller than) that of the N field, Nq. 

ii) There is only one natural frequency of oscillation, given by the masses of the fields at the 
global minimum, 

rhff, = rfiN = V2k,Nq . (22) 

This is quite different to the situation in most of the parameter regime of the standard (non- 
sup ersymmetric) hybrid potential, and it is due to the exact relation A = imposed by 
SUSY. 



in) In the limit m,^ -C k0c, and due to the fact that both fields have the same mass at the 
global minimum, there is a particular solution of the classical equations such that, 

dVo dVo 

a0"5iV' ^23) 

that is, the fields' trajectory towards the minimum of the potential describes a straight line: 

iV = ±v^(0,-0), (24) 

in the (p — N plane. That the fields follow such a trajectory or not will depend on the initial 
conditions, otherwise fixed by the inflationary dynamics and given by: 

(pin = (pc± — , Nin = — , 

/TT ZTT 

• _ ^9\^(, ^ . ^ _ , AT N 

o rr \Vini ■'■^ in) i in o rr o a r v*r in) 'in J • 

3iio 0(p 3Hq oN 

Since the initial velocities, (piniuai and Niniuai, are very small slow roll values, the fields will 
tend to move following the gradient field of the potential towards the minimum, and will end 
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in the straight hne trajectory almost exactly. In effect, we will have one oscillating mode and 
one stationary mode. The smaller the mass parameter with respect to m^, the better 
the straight line approximation. 

iv) If we take 0c Mp, we then have the relation, rh^/Ho = rfiN/Ho ~ Mp/Nq ^ 1, 
with Hq being the value of the Hubble constant at the end of inflation. This allows for 
many oscillations of each field in a Hubble time, thus the motion is not very suppressed 
by the expansion of the Universe and we should expect large amplitudes of oscillations. 
Consequently the large oscillations for a reasonable length of time will allow for the possibility 
of parametric resonance and subsequent particle production which must be investigated. 

In Ref. 1^ they have studied standard hybrid inflation with X = and found chaotic 
behaviour of the classical fields initially, that eventually becomes more regular. However, 
the SUGRA hybrid inflation model would correspond to A = where the chaotic motion 
is not present and the trajectory is given by a straight line. This can be seen in Fig. (|l]), 
where we have plotted the oscillations of the classical fields and the trajectory in the N — cj) 
plane, for the choice of parameters: 0c = 10^^ GeV , k = 10~^ and = 3.7 x 10^ GeV, which 
correspond to the choice t] ~ 0.01. We have also included for comparison the standard hybrid 
model with = g = 10~^, and the same values of 0c and m^. In both situations, the results 
will remain qualitatively the same when changing the value of 0c, and accordingly those of 
K, and m^; the main effect of reducing (increasing) the scale 0c being a weaker (stronger) 
effect of the expansion of the Universe during the first oscillations. 

Using the relation Eq. (|2^) , we can effectively replace the two oscillating fields problem by 
a single field, for example N{t), oscillating along the potential, 

V{N) = ^K^N^ + ^K^N"" - k^N^Nq. (25) 

Near the origin the potential is almost fiat, with the curvature becoming first negative and 
then changing sign and increasing as we approach and pass the minimum. Therefore, the 
field will tend to spend most of its time near the initial values, mainly varying when it 
approaches the minimum and the evolution speeds up. In a non-expanding Universe, and 
taking as initial conditions the asymptotic limit A^(— oo) = A^(— oo) = 0, the exact solution 
of the classical evolution equation is given by: 

N{t) _ 4/3 

~ 1 + 2m2t2/3 

where t = is defined when N reaches the peak of the oscillation A^(0)/A''o = 4/3. In fact the 
motion can be decomposed into two different parts, which we label region I/II respectively: 

Region I: when \N/No — 1| < 1/3 the field is near the minimum Nq, the motion is more 
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4 8 12 16 20 -1 -0.5 0.5 1 

N' (|)% 

Figure 1: Oscillations of the classical fields (left) and the trajectory in the N — (p plane 
(right). The top panels correspond to the supersymmetric hybrid model with 0c = 10''^^ GeV, 
K = 10~^ and m<^ = 3.7 x 10^ GeV. The time scale is given in terms of the approximate 
number of oscillations of the fields, N' = fh(f,t/27T, where A^' = has been defined for 
convenience as the point where the fields reach the first peak. The lower panels show for 
comparison the situation in a non-supersymmetric hybrid model with g — y/X — 10~^, and 
(pd— ^) s-nd as before. 
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harmonic and the sine squared term in the right hand side of Eq. (p6D can be neglected, i.e., 




Region II: when < N/Nq < 2/3 the field increases/decreases as N{t) oc 

In an expanding Universe, the main effect of the friction term in the evolution equation will 
be to reduce the amplitude of the successive oscillations (and thus the period), such that the 
amplitude will now decrease as the inverse of time. Then, we can approximate the behaviour 
of the N field around each peak of the oscillations (region I) by. 




where now $(t) oc 1/t. The actual decrease in amplitude over an oscillation period will 
depend on the ratio H/rh^. However, during the first stages of the oscillations the dominant 
effect will be to render the total oscillation more harmonic (symmetric) around the minimum, 
more than to reduce the factor $(t). 



3.2 Particle production: Scj), 6N. 



Due to the behaviour of the classical fields, we may expect particle production through 
parametric resonance when the field is oscillating around the peaks, that is, in region I when 
N{t)/No > 2/3. Therefore, in that interval of time we can use the approximation given in 
Eq. (1^) in order to convert the evolution equations for the quantum fiuctuations, Eq. (|^), 
into a Mathieu-like equation. Here we will be mainly concerned about the production of (p 
and particles, which indeed will occur in any SUSY model of hybrid infiation. Production 
of other species of particle will be model dependent, depending on how they couple to the 
classical fields. 



We first analyse the equations neglecting the mixing terms, Ai^f^ Eq. (^). As before, 
we will use the fact that (p and follow a straight line trajectory, and work only with the 
classical field N{t). The effective masses for the quantum fluctuations 6(j) and 6N are then 
given by: 

Ml^ = 2k'N\ (29) 
MIj, = 2k^N{2N-No), (30) 
which, when the classical flelds are near the minimum, become: 

Ml, . <{lAos[\Fj^]^,dn, (31) 




Ml^ ^ h + $(t) cos ^ m,dA . (32) 
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We have kept only the contribution hnear in the cosine and neglected the sub-dominant 
cosine squared contributions. Substituting these expressions into Eq. (|^, and defining: 



2- = / \/! (^] m,dt' , (33) 




/o V 2 V No ^ 

it is straightforward to get the parameters and qa entering in the Mathieu equation, that 



IS, 



8 f No \\, , 8 / iVo \^ qJt) 

We notice that these parameters are functions of time, as is the case in models studied in 
an expanding Universe. The time dependence purely due to the expansion appears in the 
scale factor, a{t), such that the momentum k is red-shifted as the system evolves, and in the 
damping of the amplitude of the classical oscillations, $(t), which will cause q{t) to decrease 
appreciably in a Hubble time. On top of this, we have an extra time dependence related 
to the nature of the oscillations, which enters in the factor No/N{t), and which would be 
present even in Minkowski space. 

Although the Mathieu parameters and depend on time we can infer the behaviour of 
the solutions using the stability/instability bands analysis of the canonical Mathieu equation, 
(|T2]). For a given mode k, instead of looking to a given point in the chart, we will move 
along some trajectory in the plane {Aa(t), qaif)) 0, [131' crossing different stability/instability 



regions as the fields oscillate. As long as we know the trajectory followed in the Mathieu 
chart, the standard analysis can be applied. For example, the parameters q{t) are bounded 
by, 

2$(t)> q^{t) >$(t)/2, (36) 
3$(t)> qN{t) >\m. 

so that in each oscillation we will move in the Mathieu plot downward from right to left 
during the first half period, and back again. This can be seen in Fig. (^. 

We first notice that for small values of the momentum k and $(t) ~ 1, we have that 
Afj, ~ 3g0 but v4jv ~ 2gjv The line A = 2q roughly separates the broad (to the right of this 
line in the Mathieu plot) from the narrow (to the left) resonance regimes, and accordingly we 
expect much less production of particles than particles. To see the difference between 
broad/narrow regimes, we can follow Ref. [^] in order to get the maximum allowed growth 
index in a trajectory of the type A{k) = d{k) + bq, with k being the momentum. This is 
given by, 

i In (Vl + e— + e-"^) , (37) 



max 
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Figure 2: The trajectories through the Mathieu stabihty/instabihty chart for k = showing 
how the efficiency of parametric resonance varies over each stage in the classical field oscil- 
lation. The contours are incremented in levels of /i ~ 0.1, ranging from the darkly shaded 
(stable) regions to the light (extremely unstable) regions. Modes k > will lie vertically 
above the zero mode trajectories shown. Particle production occurs mostly in the second 
instability band. 
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Figure 3: The growth index fi{A, q) for the static Mathieu Equation. 

with:, 

c=^W^. (38) 

The value /i*""^ should be understood as an upper bound on the effective growth index along 
the given trajectory. Along the line A = 2q (i.e. A; = 0) we have a constant maximum possible 
growth index, ^^"-^ ^ 0.28 p|, Q. Because we are crossing several stability /instability bands 
in only half a period, it has been argued than the effective /i will somehow be less than half 
this value, with the exact value depending on the actual time dependence of the parameters, 
and how much time is spent in each band. To the left of this trajectory, A > 2q, we have 
that c grows with q and then ^"^"-^ decreases as we move from left to right in the Mathieu 
plot. In the opposite case, A < 2q, we have the reverse behaviour with jj,"^"-^ increasing with 
large q, and we are in the broad resonance regime. 

Then, it is clear that N production is in the border line between the broad and narrow 
resonance, and the resonance can be quite efficient in producing them. With respect to 
production, we only get a non negligible value of n^""^ ~ 0.1 in a small part of the trajectory, 
as can be seen in Figs, (j^) and (^. Given the constraints in Eq. (^), at the beginning 
production occurs in the second resonance band, whilst may partially enter the first 
resonance band (see Fig. (^)). Moreover, even if the motion is not very suppressed by 
the Hubble expansion, parametric resonance takes place only during a fraction of the total 
period of the motion, so the final effective fik parameter averaging over several oscillations 
will be much smaller. Hence, for the particles we expect a narrow resonance if at all for 
momenta k < ffi^. As the damping factor $(t) gets reduced in time, A^ production will be 
shifted further to the left in the Mathieu chart, into the narrow resonance regime. 

However, A^ particles with low momenta may be produced when the classical fields are not 
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near the minimum, although this has nothing to do with parametric resonance. Whenever 
N{t)/No < 1/2 the effective squared mass of the field, Eq. (0), becomes negative and we 
enter the spinodal instabihty region. In this region, for the modes for which fc^/a^ + A^^^y < 
0, we will again have an exponential increasing solution that can be interpreted as production 
of particles. At the end of the day, it seems that production of N particles is quite efficient 
for modes with k <m^, and extremely effective for small modes below roughly 0.35m^. 

So far, so good. These would be the qualitative results if the mixing, A^^^v; "were negligible, 
but this is never the case in SUSY hybrid inflation. The analysis with the mixing term for 
the 6(f), 6N system is greatly simplified given the fact that the background fields follow a 
straight line trajectory, implying only one oscillating mode. Hence, the coupled system for 
their quantum fluctuations can also be decomposed into two orthogonal modes. 



cos 6*50 + sin 96 N , 
- sin 666 + cos 66 N . 



with tan^^ = —\/2. The time dependent frequencies for the modes 6(l)± are now given by: 

Ml = 2K^N{t){3N{t) -2No) , (39) 
Ml = 2K^N{t)No, (40) 

and the corresponding parameters in the Mathieu equation are: 

^^t) = ^^iTTi^] m, A^{t) = ^{T7i^{7-^ + ^^-^, (41) 





The Mathieu parameters for the mode appear well to the left of the line A = 2q, i.e. 
deep in the narrow resonance regime (see Fig. (|^). However, for the mode 50+ we have 
the reverse situation, and the amplification of this mode will occur in the efficient broad 
resonance regime. Moreover, the effective mass M\_ becomes negative when N{t)/No < 2/3, 
and again we will also have production of particles in region II of the classical oscillation. 



Finally, given that tan 6* ^ 0(1), both 50 and 6N are dominated by the production in the 
50+ mode, 

50 ~ - sin 6*50+ , 5A^ ~ cos 6*50+ , (43) 

leading to roughly the same rate of production for both particles. Moreover, in this case the 
mixing enhances the rate of production of N particles compared with what was expected 
without it. This enhancement is a two- fold effect: on one hand parametric resonance is more 
effective due to the mixing, on the other hand we have a greater negative frequency in region 
II, with MWmin = ~?^^/3, which will imply a wider range of momenta, k/a < 0.58m<^, in 
which particles can be more efficiently produced. 

To illustrate these points, we have numerically integrated the equations with k = 10^^, 
0c = 10^^ GeV and = 3.7x 10^ GeV, and obtained the occupation number of particles, Eq. 
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Figure 4: Initial rates of production of (j) and particles for the model with k = 10~^, 
(pc = 10^^ GeV and = 3.7 x lO^GeV. In the left side plot the results without taking 
into account the mixing term in the effective squared mass matrix are displayed; we have 
also included the oscillation of the classical field N(t) (dashed line) scaled by a factor of 5. 
The comoving momentum k is given in both plots in units of the natural frequency rh^ and 
it corresponds to the physical momentum at A^' ~ —17. In the right side plot, by the time 
particles are produced the initial value of the physical momentum k/a = 0.7m^ has been 
red-shifted below O.SSm,^. 

(|TUp. The results are plotted in Fig. and we have included for comparison those obtained 
when neglecting the mixing (left plot). Without the mixing, particles are hardly produced, 
whilst production of particles with low momenta proceeds quite quickly. Moreover, we 
can distinguish clearly in the plots the two different mechanisms of particle production. The 
peaks in Inn^ are due to parametric resonance in an expanding Universe, whereas between 
peaks Inrifc grows almost linearly in the first few oscillations reflecting the presence of a 
negative squared frequency in the evolution equation for these modes. When the Hubble 
expansion has reduced the amplitude of the oscillation such that N{t)/NQ > 1/2 (without 
mixing) or further to N{t)/NQ > 2/3 (with mixing), the occupation numbers become constant 
between the peaks. When the mixing is included (right plot), the occupation numbers for 
both particles become comparable, and the rate of production of particles is enhanced 
compared to the situation without the mixing. 



3.3 Comment on backreaction. 

Up to now, the analysis of particle production in the SUSY hybrid model has been done 
neglecting backreaction effects. Even if no other particles couple to the inflaton and A^ fields. 
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it can be inferred from Fig. that the production of their own quanta will soon affect the 
evolution of the system. Due to the presence of the tachyonic mass immediately after the 
end of inflation, these fluctuations may grow and modify the classical evolution even before 
the fields reach the minimum for the first time. After the inflaton has crossed the critical 
value, the classical fields may still move in the slow-rolling regime for a while, and we can 
consider that the Universe continues inflating for a further fraction of an e-fold Ne ~ Ht 
(for the parameters of Fig. (0) N^. ^ 0.3). The first effect of the instabilities during this 
period would be to speed up the evolution and end the slow roll regime and thus inflation 
altogether. Their contributions (50a) to the classical equations can be computed using the 
Hartree Fock approximation^ and can be interpreted as corrections to the effective masses 
of the fields. Thus, even in the slow roll regime with N{t) 0, the mass of the inflaton 
becomes, 

ml^ml + 2K'^{5N^) . (44) 

When the fluctuations grow such that both terms become comparable, the inflaton (and by 
extension the field) is pushed away from the slow roll trajectory. These contributions are 
dominated by the modes with the lower values of the momentum k, but their evolution will 
tend to follow that of the background fields (for small enough values of k, the equations of 
motion are the same). As the classical fields reach the minimum and begin to oscillate, pro- 
duction of particles in these low momentum modes practically stops, and production begins 
in slightly higher modes. This effect is similar to the "zero mode assembly" described in 
Ref. 1^ in the context of new inflation type models. Despite the differences, new inflation 
and hybrid inflation share the fact that both end with a phase transition. Under certain 
conditions, it is shown how the spinodal instabilities drive the growth of the quantum fluctu- 
ations such that they dominate completely the inflationary dynamics ending inflation. The 
super-horizon modes together with the inflaton follow the same evolution and assemble into 
an effective zero-mode field which drives the dynamics of the system. In hybrid inflation, 
the contribution due to the particles produced in the lower modes during the pre-oscillatory, 
quasi slow-rolling regime, speeds up the initial evolution towards the minimum, but it does 
not reach the level of modifying the oscillatory regime. For this to happen, we have to wait 
until higher modes have been populated and become dominant. 

On the other hand, other particles can be present and these might be more efficiently pro- 
duced during the oscillatory phase, such that they would dominate the backreaction effects. 
Given that in the original SUSY superpotential we have complex fields, their imaginary 
components, CP-odd singlet scalars, are among the particles that should also be included 
in the game. However, because the couplings to their real partners is model dependent, we 
will not attempt to study the production of these particles in a general way here. We will 
rather focus in the next section in a particular supersymmetric model of inflation motivated 
by particle physics considerations, where these questions can be addressed. 



brief description of this approximation is given later in section 4.3 and Appendix B 
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4 A particular model: (/)NMSSM Hybrid Inflation 



The 0NMSSM [1^ is a simple modified version of a well known particle physics theory, the 
next-to-minimal supersymmetric model (NMSSM). As such it is a realistic working super- 
symmetric model which supports inflation. The model has been formulated in the context 
of a SUGRA theory in |]16|. We will briefly review the global SUSY model here. 



The model is based on the superpotential: 

M/ = ANHiH2-/€N2<I' + ..., (45) 

where Hi, H2 are the standard Higgs doublet superfields, and N and $ are singlet superfields. 
As in the usual NMSSM, k and A are two dimensionless coupling constants. Since the singlet, 
$ is linear in the superpotential, the potential will be very flat in its direction enabling it to 
act as the slowly rolling inflation field. 

This superpotential admits a global U{1)pq Peccei-Quinn symmetry which forbids further 
terms and solves the strong-CP problem. The symmetry is broken at the scale of the singlet 
VEVs giving rise to an invisible DFSZ axion of that scale. Whilst the axions may provide 
a Dark Matter candidate, they must be tightly constrained in order not to over close the 
Universe or destroy the results of nucleosynthesis. Axions may be produced directly during 
the explosive preheating phase or they may be a decay product of the other massive particles 
produced during preheating. 

Previous to the start of inflation it is assumed that the field is driven to zero. During 
inflation the standard Higgs fields may be ignored and the scalar potential in terms of the 
canonically normalised real fields and N, becomes: 

^(0, N) = ViO) + ^k'N' + k'{^ - 0,_)(0 - 0,+)Ar2 + ^my, (46) 

The constant, V{0) has been added to ensure the necessary vanishing of the potential at its 
global minimum which corresponds to our vacuum and the critical values of for which the 
effective A^ mass vanishes have been written explicitly in the potential and are given below 
in terms of the soft SUSY parameters, which appear in the soft SUSY-breaking part of the 
potential, Vsoft = lm%N^ - ^kA«,0A^^ 



5c± 



1 ± 



2V2k [ \ Al 
The first derivatives of the potential are. 



1-^1- (47) 



dVo/dcP = 2«:2(0-0o)iV2 + m20, 

dVo/dN = fi:2iv[iV2 + 2(0-0e-)(0-0c+)]. (48) 
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Taking the fields to be positive, the (tree level) global minimum is found to be, 



+ 



Na 



A 



V2 



2k\ 



1 - 



4m^ 



(49) 



and thus the vacuum energy which must be zero at the global minimum is given by, 

V{0)^-Voi<Po,No) = ^^'N^. 



(50) 



In the SUGRA model |jT6[ V{0) is provided in the hidden sector by the moduli and dilaton 
contributions. 



As usual, during inflation the N field is trapped in its false vacuum, = 0, whilst the field 
slowly rolls along the potential either from above 0^+ towards the origin (hybrid) or from 
below 0c- and away from the origin (inverted hybrid |1^), depending on whether the small 
infiaton mass squared, m^, is positive or negative respectively. Once at the critical value, 
the A^ mass changes sign and both singlets roll down towards the global minimum signaling 
the end of inflation. 



The distinctive feature of having two critical values, 0c±, relies on the presence of the trilinear 
parameter A^ in the potential. Since we expect AA^o ~ 10^ GeV in order to motivate the n 
parameter of the MSSM, the inclusion in the potential of the SUSY-breaking soft parameters, 
masses and trilinears, is mandatory in our model. However, once the potential is written as 
in Eq. (^6|), it can be seen as a generalisation of the SUSY hybrid potential, independently of 
the origin of the mass parameters. In the limit Al -C |m^|,(m^ < 0), and thus 0c- = — 0c+, 
we recover indeed the standard SUSY hybrid potential, 

1-0(0, N) = ^-k\N' - 2<Plf + K^N^^ + \ml<p\ (51) 

Therefore, the main difference between Eqs. ( ^ ) and ( ^^ lies in the origin of the scale 0c. 
In standard SUSY F-term hybrid inflation, this is given in terms of an intermediate mass 
scale A coming directly from the superpotential, Eq. ([ISD , whereas in our model no explicit 
mass scale is included in the superpotential and the critical value, 0c ~ Ak/t^, is generated 
via soft SUSY-breaking. From the cosmological point of view, this in turn will mark the 
difference between the models through the value of the Hubble constant and the vacuum 
energy during inflation. 



Whatever the origin of the critical scale 0c, the parameters in the model must satisfy the 
COBE constraints and the slow rolling conditions, Eqs. (16-19). Since in our model 0c ~ 
00 ~ A"o must not exceed the bound on the axion scale, and moreover AA^o ~ ^^00 ~ 1 TeV , 
rough order of magnitude estimates for the parameters in the potential are found to be: 
(po No ^ (pc ^ lO^^GeV, A ~ K ~ 10~^°, hence the vacuum energy and Hubble constant 
are small, V{OY^^ ~ lO^GeV^, Hq ~ lO^^GeV. Note that terms involving the Higgs doublets 
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in the potential, (see Appendix A, Eq. (|65|)), will be at least a factor below the value 

of ^(0) above so essentially these fields will not contribute at all towards the vacuum energy. 
We result with a scale independent spectrum, n = 1 predicted to an accuracy of 10~^^. The 
inflaton mass, ~ eV comes from radiative corrections to its zero tree-level mass. The 
smallness of A and k, the origin of V^(0) and the massless tree-level inflaton can be explained 
in terms of a no-scale SUGRA, superstring inspired theory in where it is also shown that 



the 77-problem common to most (F-term) inflation models can be avoided here. 

In the next subsections we give the results of our investigation of the post-inflationary stage 
of the 0NMSSM model described in the previous section. Firstly we analyse the behaviour of 
the classical homogeneous background fields according to the classical potential, neglecting 
particle production (Figs. IJ-^); secondly we study the quantum fluctuations on the classical 
background of all fields in the theory, relating the growth of these fluctuations to particle 
production (Figs. |^-|T0|); and lastly we find an estimate of the point at which backreaction of 



created particles becomes important and show the immediate effect this has on the classical 
field evolution (Figs. 0-0). 



Within the approximate order of magnitude results necessary to realise the model, we have 
chosen the set of parameters below for our numerical study: 



No = lO^^GeV, K = 10-^°, = lO'^GeV, (52) 
00 V2 



More explicit details and calculations involving the other fields in the potential needed in 
the analysis of the post-inflation era are given in the Appendix at the end of the paper. 



4.1 Classical results - neglecting particle production 

Being no more than a particular version of the SUSY F-term hybrid model, the 0NMSSM 
model leads to an evolution of the classical homogeneous fields, N{t) and 0(t) (ignoring 
particle production) which follows the general pattern given in section In this case the 
inflaton field will oscillate around its non-zero VEV, 0o- The two ratios of scales that fix 
the behaviour of the oscillations are: 

(a) m^/{K(j)c) ~ 10~^^; the motion begins with extremely small slow-rolling values, and 
moreover the contribution of the tiny inflaton mass becomes negligible very soon in the 
evolution. As a result, the fields will follow the straight line trajectory (Fig. (|^) to a very 
good degree of accuracy, and the actual oscillations of the fields are indeed in phase. 

(b) Ho/{K,(f)c) ~ 10~^; the oscillations will begin with very large, lightly damped amplitudes, 
and oscillate ~ 10^ times in a Hubble time. 
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Fig. @ shows the evolution of the fields from the moment the critical point is passed. The 
time scale displayed is in terms of the approximate number of oscillations of the fields, given 
by: 

rh^t _ rriNt 

~ 2n ~ 27r ' ^^^^ 



where one e-fold {Nf. ~ Ht) corresponds to roughly ~ 2.7 x 10^ oscillation times. This 
approximation becomes more exact as the oscillations approach simple harmonic motion at 
late times in the evolution when the amplitudes are very small and the masses given by (|69|) 
become a better approximation to the actual frequency of these late oscillations. In several 
of the plots we have used the variable N' = N — N start, where A^' = is defined to be the 
moment when the fields reach the global minimum for the first time. 

We can see in Fig. ^ that the fields take a reasonable length of time, Ngtart ~ 2.9 x 10^ 
oscillation times (about 1/10 of a Hubble time) before they begin to oscillate, so in fact 
infiation continues for just a fraction of an e-fold after the critical point is passed. The fields 
then begin to oscillate with amplitudes decreasing due to the expansion of the Universe. In 
terms of a Hubble time, the total amplitude will decrease rapidly at first and more slowly as 
time progresses until after about 2 e-folds the oscillations have decreased to ~ 1/lOth their 
original amplitude. The frequency of oscillations increases as they progress, and the motion 
eventually evolves towards simple harmonic. The period of oscillation rapidly decreases 
to begin with, indicating that in comparison to the Hubble time, simple harmonic motion 
becomes a good approximation relatively early on in the stages of oscillation. However, 
particle production will occur during the first few oscillations which are far from harmonic. 

Due to the low value of the Hubble constant relative to fh^ and hence the little effect it 
has on the amplitude of the first oscillations, in the top right panel in Fig. we can 
clearly distinguish the two different parts of each oscillation: the peak, corresponding to 
more harmonic-like motion around the minimum, Eq. (pS]), and the non-oscillating section 
where the fields decay as ~ 1/t^ as they approach the initial values. This second part of the 
motion is rapidly affected by the expansion of the Universe as it causes a slight decrease in 
the amphtude, and at later stages (top right) it is further suppressed. 

The behaviour of the Hubble constant is shown in Fig. |^ where we have also plotted the 
quantity e = —H/H"^. e = corresponds to the vacuum dominated era of inflation, e = 
3/2 to the reheating phase when matter (speciflcally the decay products before they have 
thermalised) dominates, and e = 2 to the flnal (after the reheat temperature has been 
reached) radiation dominated Universe. During the oscillations, e oscillates between the 
values and 3 and its average value per oscillation, e shows the transition between vacuum 
and matter dominated states. 
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Figure 5: Classical trajectory in field space (0, A^) after inflation, (for the parameters given 
in (|52D) in the case of both hybrid and inverted hybrid inflation. The fields begin at = 
(1 ± l/v^)0o,^ = 0, and oscillate around (p = (po, N = Nq, where, neglecting particle 
production, they eventually settle. 

4.2 Quantum effects - Preheating 

Treating 6(j)a as fluctuations in the presence of the background fields, (f)a{t) which have the 
independent dynamics described in the previous section, we solve the equations of motion 
for the quantum fluctuations of all scalar fields in the theory, 

i^ka + -^^ka + MlpiPkp = 0, (54) 

where ipk = a^^'^Scpk- Here, the term e in Eq. (|^) is negligible since the Hubble constant 
is so small relative to fh^ and e ~ 0(1) (see Fig. ([^) so H is also small in this model. 
In fact, particle production will take place in a tiny fraction of a Hubble time, and we 
can safely neglect any effect due to the expansion of the Universe in the following. Since 
we are not considering the quantum fluctuations of the standard Higgs fields, 5 Hi, (see 
Appendix A), the system above reduces to two sets of two coupled equations describing the 
quantum fluctuations of the A^, 0, the axion and X fields. The axion and the X field are the 
massless/massive pseudoscalars in the model, defined in terms of the imaginary components 
of the complex fields, 

a = Nj cos tta — (pi sin , 

X = Nj sin aa + (pi cos , (55) 

with tana^ = 2(Pq/Nq. The mass squared matrices, A^^at, A^^x? given in Eqs. ( |68D and 
( [73D in Appendix A. 
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Figure 6: Top: Oscillations of the classical fields in the hybrid case at the first stage and at 
a later stage in their evolution. Bottom left: Evolution of the fields and Hubble constant 
after inflation in the hybrid case over ~ 2 e-folds, showing the damping of the amplitudes of 
oscillations caused by the decreasing Hubble constant. The fields take ~ 0.1 e-folds before 
they begin to oscillate in the shaded regions of the plot. Bottom right: The decreasing 
period of oscillations in units of (period of simple harmonic oscillations). 
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Figure 7: Behaviour of Hubble constant. Left: Small scale behaviour, showing oscillations 
in e. H/Hq is shown over the range [0.97814,0.97808], but has been rescaled to make it 
visible in the plot. Right: Large scale behaviour - Universe approaching matter domination 
(e-3/2). 
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Figure 8: Explosive production of particles, displayed here by the exponential increase in 
the occupation numbers, Uk of axions and massive X particles for k = 0.3rh^ (left); and of 
(j) and particles for k = 0.6m^ (right) with time (in units of the fundamental period of 
oscillation of the classical fields, N). The results without including the mixing terms are also 
shown for comparison. Without mixing there is no production of X particles and particle 
production is negligible. Including the mixing terms results in virtually identical numbers 
of and particles; and of axions and X particles being produced. 



The results obtained by solving the equations numerically for the parameters given in (^) 
are summarised in Figs. (§-0) where the occupation number Uk is defined in Eq. (^0]). In 
fact particle production appears to be extremely efficient. 

We have also plotted for comparison the results obtained treating each particle separately 
and solving Eqs. (|^) as if they were four uncoupled equations, ie. without the mixing 
terms (a 7^ /3) in the mass squared matrix A^^^. In particular in Fig. (H) we can see that 
with mixing there are less axions produced, but now Xs are also produced at the same rate, 
whereas N particle production increases considerably and now particles are also produced 
at the same rate. The numerical results for the (p — N system fit quite well with the analytical 



behaviour obtained in section 3.2 using the Mathieu equation 



In the axion-X system, the analysis based on the Mathieu equation will depend on the details 
of the model, for instance the value of tana^, and some of the approximations may not be 
valid, such as neglecting terms proportional to the cosine squared in the effective masses. 
Bearing this in mind, we still can get some qualitative features studying the equations 
without the mixing term. The Mathieu parameters are then given by: 



la — - 1 cos^ ctfl — sin 2tta/ V2\ ( ) ^ 0.16 I " ) , (56) 



^V-0 16 f^V 



3 \N{t) J (a2m2) 
qx ~ ^|2 + sin^Oa - v^sin^a^tana^l f^^^J ~ 0.24 



(58) 
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where the RHS numerical factors correspond to tana^ = 2. We notice that production of 
X particles would be practically non-existent, whilst axions may be produced but with qa < 
0.36, which leaves us in the narrow resonance regime. However, these are only indications 
of what could happen if the axion and X were not coupled. Due to the mixing term, we get 
that both particles are produced at the same rate, but in this case the production of axions 
will be further suppressed. 

Fig. plots an estimate of the initial rate of growth of the occupation number of each 
particle as a function of the wavenumber k, fika ~ {4:'7i)~^d(\iank)a/dN. Without taking 
into account the mixing terms, the largest rates correspond to the lowest mode, whilst with 
the mixing the effective growth index for (p and particles is negligible for k -C m^, but 
it has a maximum around k ^ 0.15m<^. As mentioned in section |3.3| , the lower modes 
follow the evolution of the homogeneous fields and their production is suppressed during 
the oscillations. For a particular typical value of k, Fig. ( pUD shows the relative rates of 
production of particles. For the beginning oscillations, the classical fields spend less than 
1/6 of the total period of the oscillation around the minimum, where parametric resonance 
takes place. In Fig. (|l3), this would correspond to the spikes in the N rate of production 
with no mixing. Thus, the relatively large values of fi{k) for low values of the momentum 
k < O.Sm^ is produced by a simpler mechanism, and is due to the presence of negative 
squared frequencies in the evolution equation of the quantum fluctuations, as can be seen 
in Fig. (0). On the other hand, we want to stress the fact that with mixing the rate of 
production for the axion and the X are brought together, but this now corresponds to a 
smaller rate of production for the axion. We will see in the next section that indeed the 
total number of axions produced is several orders of magnitude smaller than that of the 
infiaton and N particles. 

The situation for the pseudoscalars may change depending on the model. For example, in 
models where 0o = 0, (0c+ = —4>c-) there would be no mixing among the massive and 
the massless pseudoscalars and the axion would be the imaginary component of the N field. 
When particle production begins, the effective squared mass for the axion would be given 
by (see Eq. (|70[) in Appendix A), 

Ml = 2k^N{N - No) . (60) 

which again becomes negative whenever N{t) < No- In this case, we would find that indeed 
the production of axions may proceed at the same rate as that of the infiaton and A^ particles 
in a extremely efficient way. 



4.3 Effects of backreaction. 

Due to the extremely large numbers of and A^ particles produced during the first few 
oscillations of the background fields, these particles will quickly give a non-negligible contri- 
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Figure 9: Typical values of the growth parameter, /i^ over the first few {N' < 40) oscillations 
of the classical fields as a function of the momentum k in units of m^. 




Figure 10: Comparison of initial rates of production of all particles for a typical momentum, 
k = 0.3m^. 
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bution to the evolution equations. Therefore, instead of the tree-level equations considered 
until now, the full dynamics taking into account the interactions among the particles being 
produced should be included in the evolution equations. In practice this means that the 
polarization tensors for the different fields must be computed beyond the tree-level approx- 
imation. A consistent treatment of all non-linear effects involved can be obtained by using 
lattice techniques R B]. However, in order to get some insight on how backreaction affects 



the process of preheating, the Hartree approximation pO| , |2T|] for the potential should be 
good enough. 



Within the Hartree approximation, the main one-loop contributions to the tadpole and mass 
terms in the potential are included, not only for the background fields but for any quantum 
field. A summary of how to expand the potential and the new expressions for the effective 
masses are given in Appendix B. Including one-loop corrections, the effective mass of all the 
fields pick up new non-thermal contributions proportional to, 

= / ^i^0^«r- (61) 

For example, the effective mass for the inflaton field (j) and its quantum fluctuation will 
become: 

Ml^ = 2k\N^ + {6N^) + cos^ aa{{6a^) + t^{6X^) + 2t{6a6X)) , (62) 

where t = tana^ = 2(I)q/Nq is the mixing angle in the axion-X sector defined in Appendix 
A. The average can be seen as a closed loop of a particle (50q,. The Hartree approxi- 

mation does not cover all possible types of contributions to the polarization tensors, and the 
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effects of rescattering among the different species are missed. However, by integrating the 
coupled system of equations together with Eq. (0), it succeeds in summing up the infinite 
set of diagrams of the "cactus" or "daisy" type. The overall effect of the change in the 
effective masses will be to shut down the explosive production of particles. On one hand, 
it will increase the frequency of the oscillation of the background fields, hence decreasing 
their amplitudes; on the other hand, it will render the squared frequencies of the quantum 
fluctuations positive, ending the process. 

Using the results of the previous section we can estimate the time at which backreaction 
becomes important, that is, when (50^) ~ 0(00) ~ 0{N^). The expectation values are 
related to the total number of particles produced, and the later can be estimated given the 
growth index function fi{k). It is then clear from Fig. (6) that we have to worry mainly 
about and production, and that the total number of axions and X fields is several orders 
of magnitude smaller. This simplifies life, because the function n{k) for the fields and N 
can be quite nicely approximated by: 

/ \l/3 

fi{k) ^ 0.6 — exp{-5{k/m^f/'^) , (63) 

valid for the first few oscillations. The total number of particles (or A^) produced at a time 
N' = rh^t/{27r) is then: 

n^N') = J dkk^e'^^^^^""' , (64) 

and we get ^ C'(0o) when A^' ^ 20. Therefore, backreaction effects will become important 
just after the third oscillation of the inflaton field. Soon after that, the results given in the 
previous section are modified and quantum fluctuations cease to be produced at a exponential 
rate. As a consequence, we will never reach the point where we have produced too many 
axions (see Fig. (|])). 

In Fig. d^) it can be seen that only particles with low momenta, k < m^, are effectively 
being produced. Then, the value can be used as a cut-off in the momentum integrals, 
especially when computing the expectation values (50^). In principle, these quantities are 



quadratically divergent and have to be renormalised ||2^. In practice, we avoid this subtle 
point using the cut-off; the rest of the couplings and mass parameters appearing in the 
equations are then taken to be the renormalised ones. 

The numerical results confirm our estimation. In Fig. (0) it is shown how the oscillations for 
the background fields change after A^' ~ 20. However, they are still in phase and follow the 
trajectory given in Fig. (Kf). Recently, it has been stressed that, in models with symmetry 



breaking, backreaction of the quantum fields could restore the original symmetry ||23|, if 
they come to dominate the effective mass during the stage of preheating. This is not the 
case in our particular model. Backreaction itself prevents the quantum fluctuations from 
reaching the stage where they might restore the symmetry. However, the shape of the 
potential changes and when particle production has effectively stopped, the classical fields 
end oscillating around a minimum different to 0o, Nq. We stress that the whole preheating 
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Figure 12: Oscillations of the classical fields with backreaction effects included. 



process takes place within much less than a Hubble time. This means the expansion of the 
Universe can be neglected during this time. The decrease in the amplitude of the oscillating 
fields in Fig. (|1^) is purely due to the presence of the fluctuations. 

In Fig. ( |I3| ) we give the rate of production for all the different species of particle with a 
typical momentum k = O.Srh^. Soon after A^' ~ 20 they reach a plateau, signaling the near 
end of the preheating period. Fig. ( p!^ shows the total number of particles per comoving 
volume being produced. We end with the ratio of scalars to pseudoscalars being O(10^°). 
Preheating is very efficient in producing and N quanta, and in fact it is this effect which 
precludes other particles from being so efficiently produced. Backreaction due to and N 
affects the evolution of every field (classical and quantum) in the model. 

This, however, may not be the end of the story. The Hartree approximation neglects the 
effects of rescattering among the particles produced, but it is possible that these terms could 
be the same order of magnitude as the expectation values [§]. Lattice simulations 

1^ for one-field inflation models have shown that the maximum value of {S(f)'^) is overesti- 
mated in the Hartree approximation when the resonance is very efficient. However, once 
this maximum value is reached the resonance ends qualitatively in the same way as in the 
Hartree approximation. It is difficult to extrapolate from these studies to our model. We 
may have overestimated the total number of (f) and produced, and/or it may happen that 
rescattering may induce an extra little production of axions and X. But it is unlikely this 
could reduce the large gap among and n^. We can venture that, after preheating, the 
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Figure 13: Hartree approximation: Rates of production of all particles for a typical comoving 
momentum k = O.Srh^. 

Universe will end with the classical fields still displaced from the global minimum, and fill 
up with nearly non-relativistic particles, a large fraction of them being and A^. A fraction 
of the energy storage in the inflaton field at the end of inflation has been transferred to the 
particles produced, and this will set the initial conditions for the subsequent evolution of the 
system, as given by the standard theory of reheating; in the Hartree approximation roughly 
60% of the initial energy density is converted into quantum fluctuations. Moreover, by the 
time reheating is complete, any density of stable particles (i.e. axions) would be diluted by 
the entropy released in the decay of the and N. 



5 Conclusion 

We have analysed in this paper the preheating era which occurs immediately after the end 
of inflation in a general SUSY model of hybrid inflation. From the point of view of the 
inflationary potential it is no more than a particular case of the general hybrid model, 
Eq. (p!3D. Still, preheating in this particular model is found to be quite different when 
studied in more detail. Supersymmetry leads to equal masses at the global minimum for 
the real singlets relevant to inflation (0 the inflaton and A^), and therefore to the presence 
of a common frequency of oscillation for the classical fields. That is, the supersymmetric 
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Figure 14: Hartree approximation: total number of particles produced. 

potential allows an approximate solution with only one mode of oscillation for the classical 
fields, the orthogonal one being stationary. This solution becomes exact as the explicit 
tiny mass of the infiaton, m^, tends to zero. Due to the inflationary dynamics (slow roll 
conditions) this mass is indeed always negligible, independently of the value of any other 
parameter. Moreover, the slow roll scenario also fixes the initial conditions at the beginning 
of the reheating era such that we begin with very small initial velocities. Therefore, the 
fields will follow almost exactly a straight line trajectory, oscillating with equal frequencies 
and proportional amplitudes. Hence, for the purpose of studying particle production, the 
SUSY hybrid model can be simplified to a single field model (the oscillating mode), with the 
potential given in Eq. (pSf). 

Since we are dealing with an effective single field oscillating, we can attempt to convert the 
evolution equations for the quantum fluctuations into Mathieu-like equations, whose solu- 
tions in term of stability/instability regions is well known. We have first studied production 
of and particles, which is completely independent of the details of the model. Even 
in a first approximation neglecting the effects of backreaction, the evolution equations are 
coupled through mixing terms in the effective mass matrices. The equations for S(j) and SN 
taken independently (that is, neglecting the mixing term) result in a higher rate of produc- 
tion of particles than 0, the latter being almost completely suppressed. However, the 
effect of the mixing term is to bring together the rate of production for both particles, and 
overall to become more effective indeed. The fact that we achieve almost the same number 
of particles is general to any system of coupled particles with non negligible mixing. In other 
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words, during the preheating proccess we can produce particles which in principle are not 
expected to exist, if they are coupled to other particle species which are efficiently produced. 
The rate of production of the latter will be increased/decreased depending on the strength 
of the mixing. 

Given that hybrid inflation ends by a phase transition, for some values of the fields tachyonic 
masses will be present, for example that of the field. As well as an initial production 
of particles when the fields are ffist moving towards the minimum, if the amplitude of 
the oscillations is large enough there will be additional production of particles between 
oscillations when the classical fields are in the unstable region of the potential with negative 
curvature. Therefore, it is the combination of parametric resonance and tachyonic mass 
which causes the production of and particles to be so efficient. 



Production of (f) and A^ fluctuations will rapidly give rise to non-negligible backreaction 
effects which will modify the subsequent evolution of the system. However, other particles 
may be efficiently produced in the model which can also contribute to this effect. With 
this issue being model dependent, we have proceeded to analyse the situation in a particle 
physics based model of inflation, the 0NMSSM [1^. This model is not only suitable for 
inflation, but also solves the strong CP problem by an approximate Peccei-Quinn symmetry 
and leads to an invisible axion, which could be dangerously produced during preheating. 
However, in this case the mixing term between the axion and the massive pseudoscalar 
slightly suppresses the axion production rate, which turns out to be much smaller than 
that of the real scalars. Nevertheless, we obtain the same general result of comparable 
rates of production for both pseudoscalars. Backreaction effects are then dominated by the 
production of the real components of the fields, which quite rapidly succeed in preventing 
the production of any other particle species. Therefore, axions in this model are kept safely 
under control, thoughtout and after the preheating era. 
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A Appendix 

Here we present some further calculations on the 0NMSSM model, involving all scalar fields 
in the theory. 

The superpotential, 

W = ANH1H2 - kN^^ + 
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leads us to the full Higgs part of the tree level scalar potential: 



\/o(Hi,H2,$,N) = /€2|N|^+(A2(|Hi|2+|H2n+4/€2|$|2 + ^^)|N|2 

-KA^{m^ + $*N*2) - 2kA($*N*HiH2 + h.c.) + A^lHiHa]^ 
+m?|Hip + m2^|H2p - A^a(NHiH2 + N*HtH;) (65) 
+mj|<l>p+ Fz)(Hi,H2), 

where Vd(Hi,H2) is the usual D-term of the MSSM. The full spectrum of scalar fields is 
listed below, 



$ = (0 + i(f)i) 
N = {N + iNj)/V2, 
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Hi = (;i) = ^( (66) 
H2 = 



with VEVs given by: 



($) = </>o, (N) = No, (Hi) (H2) = f ° ) , (67) 



where the electroweak values, i/i ^ 1/2 250 GeV. 



The full 12 X 12 mass squared matrix for the scalar fields is defined by, M^/s — d'^V/d4>ad4>f3, 
where (pa is taken to mean each of the real and imaginary fields in (^). When the standard 
Higgs are at their VEVs, Ai^ consists of two lots of 4 x 4 matrices for the CP-even and CP- 
odd neutral sectors, and two standard MSSM 2x2 blocks for the charged mass matrices. 
All other entries are zero since the charged sectors do not mix with the neutral ones and the 
CP-odd and -even sectors do not mix. 

In fact the 4x4 matrices can be further split into two 2x2 matrices since mixing between the 
singlets and the standard neutral Higgs (at least when taken at the VEVs) are suppressed at 
the very least by a factor 0(k) with respect to the other terms, so can be neglected. (Note 
that in general though, these terms may be non-negligible and then must be included in Eqs. 
(^, leading to the possibility of parametric production of Higgs particles and also a change 
in the production rates of other particles.) We are then left with standard MSSM matrices 
for the CP-even neutral Higgs and also for the CP-odd sector with the usual massless mode 
and mixing angle given by, tan/5 = ^2/1^1- The mass squared matrix for the real singlet fields 
is given by, 

2 _ ( 2k'N' + m% AK'N{(j> - 0o) ^ 

■^'^^-Wn{<p-<Po), 3^'N' + 2k\<P-<P,^){<P-<P,.))' ^^^^ 

where (j) and N denote the real classical fields, the imaginary parts of which are held at zero. 
Again, here terms due to coupling to the standard Higgs are very suppressed with respect 
to the other terms so have been neglected. 
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At the global minimum, ( |68D is diagonal with the masses equal, 

ml = ml = 2k^NI (69) 
where the bar is used to stress that the masses are at the global minimum. 

Similarly the CP-odd mass squared matrix for the pseudoscalar singlet fields, 0/ and Nj is: 

■^^^jNj - [ Ak^(I)oN, k^N^ + 2/^2(0 + 0,+)(0 + 0,_) )■ 

As in the CP-even matrix (|68D, we have neglected the small couplings to the standard Higgs. 
Diagonalising when the fields are at their VEVs gives us the massless mode (axion, a) rha = 0, 
and the massive pseudoscalar, X, fh\ = 2k2(40q -|- A^^q), with mixing angle, given by, 

tana, = 20o/Aro, (71) 

with the axion and it's massive counterpart defined in terms of the imaginary components, 
0/, Ni, by: 



a = Nj cos aa — (pi sin a,, 
X = A/ sin -|- 0/ cos a,. (72) 
The mass squared matrix, in the basis {a,X}, then becomes, 

Mlj, = 2 cos' aa^^ (73) 

where, 

F = Yi + Y^t" - 2Y3t, 

G=iY,-Y2)t + Ys{l-t'), 

H = Y^e + Y2 + 2Yst, 

Fi = K'Ny2 + k\4> + 0e+)(0 + 0c_), 

.2 7\7-2 I ^2 



Y2 = K'N' + m'^/2, 
Ys = 2^<P^N, 
t = tanoa. 



B Appendix 



With particle production in mind, we now proceed to obtain the equation of motion in the 
Hartree approximation. Written explicitly in terms of the four fields, 4>, N, a and X, the 



full potential, (^61) becomes: 



cos^ aaiAa'^ + Ba^X + Ca'X' + DaX^ + EX^) 

cos'' aa{Fa' + GaX + HX'^), (74) 
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where, 



A = {i/A-e), 

B = t(2t^ - 1), 
C=(t4- (5/2)^2 + 1), 

D = t{2-t^), 
E = ^2(^2/4 + 1), 

and F, G and H are as given in Appendix A above with 0, replaced now by (p, N. The 
fields are decomposed as in Eq. ([|), and the potential is expanded assuming the factorization: 

6(j)aS(f)p6(f)^6(f)ri {6(f)aS(j)/3)6(f)^6(j)rj + permutatious + Constant . (75) 

Each term in the factorization can be seen as the contribution from a one-loop diagram, 
with the contractions (■ ■ ■) accounting for the internal propagator in the loop; only those 
terms which correspond to physically allowed diagrams are included in the expansion. In 
this way, one-loop contributions to the tadpole and masses are included, such that we end 
with a potential quadratic in the quantum fluctuations. Using this potential and the tadpole 
constraint {6(j)a) = we obtain the equations of motion as given in Eq. (2) and Eq. (5). 
The first derivatives dV/d(j)a, including the additional terms, are now given by: 

(^] = ^ + 2K\<f)-<f)o){6N^) + AK^N{6N6<f)) 

+2K\(j) + 0o) cos^ aa{{6a'^) + 2t{6a6X) + t^{5X^)), (76) 
+ K^N{3{6N^) + 2(502)) + AK\(f) - (f)o){SN6(f)) 



Hartree 



\dNL. , ON 

\ / Hartree 

cos^ aa{F'{6a^) + G'{6a6X) + H'{6X^)), (77) 

with dVo/d(p and dVo/dN as given in Eqs. (|5), and we have defined, 

F' = {2t^ + 1)A^ - 4t0o, 
C = 40o(l - f^) - 2tN, 
H' = {t^ + 2)N + 4t0o- 



Eqs. ( pBD and ( [73| ) for the mass squared matrix also obtain new contributions such that: 
A^2 _ , .2 , 2 / 2{5N'^) A{5(t)5N) \ 

■^'f'NHartree " ■^<I>N + ^ [^{S^SN) 3{SN^) + 2{6(P^) ) 



,2 „„„2 



, , 2{{6a^) +2t{6a6X) +t^{6X^)) 
"^"^ ^ (2^2 + 1) (5a2) - 2t{5a5X) + {f + 2) {5X^) 



<2 ,.2 , ,.2_2„, / (l + 2t2)(<5Ar2) + 2(502) 2t{6<p^)-t{6N^) 



Max Hartree " -^aX + COS ^ ^^(^^S^ _ ^^^^2^ (2 + t^){SN^) + 2^2(502^ 
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2 4/ 12A(5a2) + QB{5a5X) + 2C{5X^) 3B{6a^) + AC{6a6X) + 3D{6X^) \ 
+K COS aa n,B(5a^) + AC{5a5X) + ?,D{5X^) 2C{5a^) + 6D{5aSX) + 12E{SX^) J " 

(79) 



In addition, the Friedmann equation for a Universe containing classical fields and particles 
now reads: 



in 



3Ml 



^0) 



with pp being the energy density of the different particles produced, i.e., 

k^dk 
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